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Nonlinear photon transport in a semiconductor waveguide- cavity system containing a 
single quantum dot: Anharmonic cavity-QED regime 

S. Hughefl and C. Roy 

Department of Physics, Queen's University, Kingston, Ontario, Canada K7L 3N6 

We present a semiconductor master equation technique to study the input / output characteristics 
of coherent photon transport in a semiconductor waveguide- cavity system containing a single quan- 
tum dot. We use this approach to investigate the effects of photon propagation and anharmonic 
cavity-QED for various dot-cavity interaction strengths, including weakly-coupled, intermediately- 
coupled, and strongly-coupled regimes. We demonstrate that for mean photon numbers much less 
than 0.1, the commonly adopted weak excitation (single quantum) approximation breaks down, even 
in the weak coupling regime. As a measure of the anharmonic multiphoton-correlations, we compute 
the Fano factor and the correlation error associated with making a semiclassical approximation. We 
also explore the role of electron-acoustic-phonon scattering and find that phonon-mediated scatter- 
ing plays a qualitatively important role on the light propagation characteristics. As an application 
of the theory, we simulate a conditional phase gate at a phonon bath temperature of 20 K in the 
strong coupling regime. 



PACS numbers: 42.50.Ct, 42.50.Pq, 78.67.Hc 



I. INTRODUCTION 

The ability to couple waveguides and cavities offers 
exciting opportunities for integrated quantum optical de- 
vices using solids^^. In particular, planar photonic crys- 
tals offer a technology platform, where quantum bits 
(qubits) can be manipulated from quantum dots (QDs) 
placed at field antinode positions within the cavity or 
waveguide^^. Integrated semiconductor micropillar sys- 
tems also show great promise for quantum optical appli- 
cations^!^, working at the few photon level. 

Recently, there have been several successful demon- 
strations of coherent light propagation effects in various 
semiconductor systems, including planar photonic crys- 
tals and micropillars. Bose et ai^ measured the exciton- 
induced doublet (polariton splitting) through waveguide 
mode transmission in a photonic crystal waveguide- cavity 
systen>iS [cf. Fig. 1(a)], while Loo et ali^ probed the 
strong coupling in a micropillar via coherent reflection 
[cf. Fig. 1(b)]; Young et alM demonstrated first steps to- 
ward a conditional phase gate using light reflection from 
a micropillar. Common to the analysis of all of these ex- 
periments has been the application of the weak excitation 
approximation (WE A) , where at most only one quantum 
is assumed. For example, Ref. IH suggested that their 
experiments were likely at the "single photon level" for 
less than 0.1 photons per cavity lifetime, so they applied 
a WEA solution. The same assumptions are tacitly made 
by many other groups, in excitation regimes where the 
mean photon number is well below those associated with 
saturating the QD exciton^. These useful formalisms 
have been very successful and certainly help to clarify 
the basic physics of low-intensity photon transport. 

For increasing field strengths, however, the validity for 
the WEA becomes questionable, and there can be quan- 
tum nonlinearities in the system due to multiphoton cor- 
relations. Giant optical nonlinearities were studied by 



Auff eves- Gar nier et alt^; their semiclassical approach 
adiabatically eliminated the cavity mode and included 
effects outside the WEA (the "Purcell regime"); natu- 
rally with such a semiclassical approach, there is no in- 
fluence from the higher lying levels of the anharmonic 
Jaynes- Gumming (JC) ladder, so it cannot be applied 
in the strong coupling regime. Most photon transport 
approaches also neglect the details of electron-acoustic 
phonon scatteringi^ — apart from the inclusion of a 
Lorentzian decay rate for the exciton, i.e., broadening of 
the the zero phonon line (ZPL). Recently, several works 
have shown that coherent excitation of semiconductor- 
QD systems can easily go into the anharmonic cavity- 
QED (quantum electrodynamics) regime^SiSl. Thus the 
questions arise: (^) To what extent can one safely employ 
the WEA for these emerging semiconductor waveguide- 
cavity and micropillar systems? (ii) When are multi- 
photon correlation effects important (quantum nonlin- 
ear regime)? (Hi) What is the role of electron- acoustic 
phonon scattering and how does this mechanism differ 
from a simple-Lorentzian pure dephasing model? 

In this work we attempt to answer these questions 
for the regime of coherent photon transport, and show 
that generally one must include multiphoton effects 
and phonon scattering within the theoretical formal- 
ism. Even for weakly-coupled QD-cavity systems, higher- 
order quantum correlations effects are shown to be signif- 
icant. Our paper is organized as follows: In Sect. [Ill we 
introduce the general theoretical technique to simulate 
coherent photon transport outside both the WEA and 
the semiclassical approximation. The theory is based on 
a quantum master equation (ME) formalism where cav- 
ity and dot decays are introduced by Lindblad superop- 
erators, and phonon interactions are included from an 
effective phonon ME^^ that includes electron-phonon in- 
teractions at a microscopic level — derived using a polaron 
transforn>2£i2^. In Sect. IIIH we first present calculations 
with no acoustic phonon coupling, for several different 
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FIG. 1: (Color online) Schematic of two semiconductor cavity- 
QED systems containing a quantum dot: (a) waveguide- 
cavity system (e.g., made up of a planar photonic crystal) 
where we assume ki = Kr = "2 kq (where kq is the out of plane 
loss), and (b) a micropillar cavity, (c) Dot and cavity energy 
levels, where the exciton system (|+/— )) interacts with an 
acoustic phonon bath (see text). 



QD-cavity couplings; both the weak coupling regime and 
the intermediate-to-strong coupling regime are investi- 
gated. We demonstrate that substantial deviations from 
the WE A can result, even for very small mean photon 
numbers <C 0.1. We then modify the ME approach to 
include the mechanism of electron-acoustic-phonon scat- 
tering, and study the impact of electron-phonon inter- 
actions on incoherent scattering and on coherent renor- 
malization of the exciton-cavity coupling rate; qualita- 
tive differences from a simple Lorentzian decay model 
are found. As an application of the theory, we study the 
transmission of light in the strong coupling regime and 
simulate a conditional (exciton-induced) phase gate. We 
conclude in Sect. IIVI 



II. THEORY 

We wish to describe light propagation for a QD-cavity 
geometry, where the input and output fields can be iden- 
tified separately from the cavity region in which the QD is 
assumed to be embedded. An example waveguide- cavity 
system is shown schematically in Fig. 1(a). For a con- 
tinuous wave (cw) waveguide mode of a photonic crystal 
system, the classical WEA reflectivity was previously de- 
rived by Hughes and Kamad a^^i^^ : 
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(1) 

Fo = 2a>:o is the cavity decay rate through vertical scatter- 
ing {unloaded cavity broadening), Fc = 2f<ic = 2(/^^ + k.^) 
is the cavity- waveguide coupling rate (which is inversely 



where the self-energy cj5](cj) = ujg'^/{uJx 



proportional to the group velocity of the waveguide 
mode^), ujc is the cavity mode resonance, uj^ is the target 
exciton resonance of the QD, and F^ = 7 + 7MS the to- 
tal decay rate of the exciton — including radiative (7) and 
non-radiative, pure dephasing (7') processes. The total 
cavity decay rate is F^ = Fq + Fc and the exciton-cavity 
coupling rate, g oc d'^/Vef^, where d is the dipole moment 
of the exciton and Vef[ is the effective mode volume. The 
corresponding transmissivity is simply t = 1 + r, and one 
can also define the reflection and transmission, through 
R = |rp and T = \t\^. Similar expressions have been de- 
rived by other groups, e.g., Refs. fT3l23l25l . With the dot 
on resonance with the cavity, then a polariton doublet 
coincides with the vacuum Rabi splitting which can be 
observed in transmission or reflection; this normal mode 
doublet can occur even if the dot is not in the strong cou- 
pling regimei^i^iiSS, though ultimately the doublet fea- 
ture is lost at high temperatures due to phonon bath 
coupling^^. The analysis of coherent reflection from a 
micropillar system is similar; one has r^pin = 1 — y^Tpc, 
where 77 is a measure of in/out coupling efficiencysi^, and 
one also makes the following replacements: Fc Tq (ver- 
tical scattering) and Fq Tg (sidewall scatter). 

The above WEA formalism does not distinguish be- 
tween radiative and pure dephasing processes of the 
QD exciton (which a semiclassical approach can); nor 
does it take into consideration the acoustic phonon bath 
or multiphoton effects. Ignoring the complexities of 
phonon scattering for nowi^i^i^^ — which we will intro- 
duce below — these analytical formulas are expected to 
work for only weak excitation conditions (strictly linear), 
and for pure dephasing rates^, 7' <C ^^F^/[4(cjc — + 

(r*/2)2]. 

To go beyond the WEA, and the semiclassical ap- 
proches, we will use a ME approach where exciton- 
photon interactions are easily included to all orders. Re- 
ferring to Fig. IHa), we relate the left/right output oper- 
ators to the cavity mode operator through^'' 




{ain{t)) + V 

'2^ (a(t)) , 



(a(t)), 



(2) 
(3) 



where, for a coherent cw input state, (ain) = 
ir]c/ (2V2^c), with r]c the cavity pump rate. Following 
the solution of the ME (discussed below), the steady-state 
transmissivity and reflectivity are obtained: t = \t\e'^^^ = 
Rut)ss/(«in) and r = |r|e^^- = (aj,^t)ss/(«in) , where 
(j)t and (j)r are the phases. Working in a frame rotating 
with respect to the laser pump frequency, ul, the model 
Hamiltonian can be written as 

H = hAxLO'~^o-~ + hAcLd^CL + hg{a~^a + aV~) 

q q 

where hq{h]^ are the annihilation and creation operators 
of the phonons, a is the cavity mode annihilation opera- 
tor, are Pauli operators of the electron- hole pair 
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(exciton), /\ciL = — = x^c) are the detun- 

ings of the exciton (ujx) and cavity {uJc) from uj^^ and 
^ drive — ^Vc{ci + CL^) dcscribcs the coherent cavity drive 
(excited through the waveguide channel). 

To obtain the ME, we first transform Eq. (|4]) to a po- 
laron frame, which formally recovers the independent bo- 
son model (IBM)^— ; the IBM is known to accurately 
describe the characteristic lineshape of a single exciton 
coupled to a bath of acoustic phonon^. Specifically, the 
polaron M E^^i^^i^-*^ introduces coherent electron-phonon 
coupling exactly, while incoherent phonon interactions 
are treated at the level of a second-order Born approx- 
imation. Further details og the model are discussed 
in Refs. [2QI131I . The time-convolutionless ME takes the 
form^^: 

|^ = l[i7;„p(t)] + + /:ph(p), (5) 

where the polaron-transformed system Hamiltonian is 
i/^y, = h{A,L - Ap)a+a- + hA^La^a + {B)Xg + H^^,^^, 
with (B) = exp[-^ duJ{co)/u^ coth{f3nuj/2)]{p = 
l/kbT), Xg = hg{a^a~ +<j+a), and Ap = duj J {uj) / uj . 
In what follows, we will absorb the polaron shift 
(Ap) into the definition of ujx- The phonon spectral 
function2£i^i^, J(^uj) = o;^ exp(-cj^/2cj^), describes 
electron-acoustic phonon interaction via a deformation 
potential coupling. 

Using a Markov approximation, the incoherent phonon 
scattering term is defined as^ 

1 / 

CM = --2 (^™(^) 

"'O m=g,u 

X X^,e-^<^"/^X^e^<^"/V(t)l +H.C.) , (6) 
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Green function^i^: Gg{t) 



where Xu = 
are the polaron 

(5)2 (cosh[(/)(t)] - 1) , Gu{t) = (5)2 sinh[0(t)], with 
0(t) = duj^ [coth{f3huj/2) cos{ujt) - i sin(cjt)]. In 
Ref. [22I, an effective Lindblad ME has been shown to 
yield very good agreement with the full polaron ME solu- 
tion above. In this way, one defines the phonon- mediated 
incoherent scattering processes through 

^ph(p) = %^^(^+a) + ^^/:(ata-), (7) 

where C{D) = 2DpD^ - D^Dp - pD^D, and the scatter- 
ing rates are obtained analytically, fron>^ 



(8) 



where Acx = ujc — ^x is the cavity-exciton detuning. The 
rate T^i^ describes the process of cavity excitation and 
the emission of a exciton, via phonon-induced scattering. 



of a cavity photon. For completeness, one can also add 
in the Stark shifts^^, but these are found to be negli- 
gible for the regimes considered in this work. An al- 
ternative weak-phonon-coupling theory is presented in 
Ref. Since this Lindblad-ME-form is considerably 
easier to work with and helps to identify the physics of 
the scattering processes in a more transparent way, we 
will use this latter form for our phonon calculations in 
this paper. 

Phenomenologically, we also include Liouvillian super- 
operator^^^i^: 

C{p) = —{2a~pa~^ — cr^cr~ p — pcr^cr~) 



(2apa^ 



ap - 



pa^a) 



— {(Jiipan - (Tnanp - panan), 



(9) 



where 2tvt = "^i^o + 2/^^ = F^, 7 = 7(5)^ is the radiative 
decay rate of the exciton, 7' is the pure dephasing rate 
of the exciton, and an = a~^a~. Figure 1(c) shows a 
schematic of the model, which we solve in a basis of \n) 
photons and the two QD states |— /+). The WEA cor- 
responds to a three state model, i.e., only including the 
ground ground state and and the first two ladder states 
of the JC model. 

For suitably large input fields, many JC ladder states 
may be involved. To better highlight the role of quantum 
statistics in the system and to help quantify the possible 
failure of a semiclassical approximation, we introduce a 
function for computing the relative correlation error in 
the expectation (aV~)gg compared to (a^)gg {a~)^^, 



CE 



l(^^^-)ss-(^nss(^-)ssl 
l(«^^-)ssl 



(10) 



where 100 CE is the percentage relative error; in a semi- 
cassical picture, of course CE = 0. We also compute the 
mean photon variance, defined through the Fano factor: 



{ia^ar),,-{aU)l 



(11) 



and ry 



describes exciton excitation via the emission 



III. RESULTS 

A. Photon transport with no acoustic phonon 
coupling 

We first investigate the system with no acoustic 
phonon coupling — apart from the ZPL broadening (i.e., 
no £ph process and (B) = 1). For calculations, we 
use material parameters that correspond closely to those 
in experiments: 7^; = 1 /ieV, tZc = "^i^o = 50 /ieV, 
7' = 4 /ieV (unless stated otherwise), and we will study 
various values of g and r]c. 

In Fig. [2] we study the weak-coupling regime, with 
g = 20 peV ^ 0.32 tzt^ using a fairly weak excitation 
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FIG. 2: (Color online) Transmission characteristics for 
a weakly- coupled QD-cavity-waveguide system, with g = 
20 /xeV ~ 0.32 ^tt and r\c — 0.5 g. (a) Transmission with 
(dashed) and without (solid) the WEA approximation. All 
other calculations in (b-d) use the full multiphoton approach, 
(b) Steady-state cavity (solid) and exciton (cavity) popula- 
tions, (c) Relative correlation error of (ttV")^^ if a semiclassi- 
cal approximation was used, CE, and the Fano factor (photon 
number variance), (d) Cavity photon moments ((tt^)^^^)^^ 
versus m. The other parameters are as follows: 7 = 1 /xeV, 
J — 4/ieV, and nc — 50 /xeV (/^t = 62.5 iieV). Of note, the 
polariton doublet in (a) is not due to strong coupling, since 
we are in the weak coupling regime for these simulations (see 
text). 
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FIG. 3: (Color online) As in Fig. [21 but with the larger pump 
field r\c = 2.5 g. We clearly see that multiphoton correlations 
are needed, even for this weakly coupled system {g <^ Kt). 



field of r]x = 0.5 ^. This field value was chosen to be 
small enough that the cavity population is significantly 
lower than 0.1, but large enough to see a breakdown of 
the WEA. We have also confirmed that this value of 
g/f<i yields no vacuum Rabi splitting. In Fig. [2fa) we 
show the transmission versus detuning with (dashed) and 
without the WEA (solid); a few comments are in order: 
(i) we confirm that the polariton doublet appears even 



thou gh w e are not in the strong coupling regime (also see 
Refs. I14II23I26I ): (ii) the WEA breaks down with already 
qualitative differences of more than 40% near ^ cjc; 
(Hi) with the chosen value of g^ the region of "trans- 
parency" is notably very weak, which is a consequence 
of the finite QD broadening (through 7 and 7'). This 
latter observation can be contrasted with the results this 
can be contrasted to the results in Ref. l14 where such 
broadenings were not included; these ZPL broadenings 
are essential to include for a realistic QD system. In 
Fig. [2fb), we show the (numerically-exact) exciton and 
cavity- mode populations, confirming that the largest cav- 
ity population is well below 0.1; however, we note that the 
fundamental condition for the WEA is not a low number 
of photons, but a negligible excitation of the dot — and 
the exciton population is evidently no longer negligible. 
In Fig. [Sfc), we display the correlation error, CE, and 
the Fano factor, F. Even at these weak drives and small 
g (weak coupling regime), it is clear that a semiclassical 
approximation can fail, especially near uj ^ Uc where the 
percentage error of assuming (aV")^^ ~ {^^)ss (^~)ss 
as much as 50%. In Fig. [2jd), we plot the corresponding 
cavity-photon moments {{a^)^a^)^^ for the first few pho- 
ton states, which suggest that only the first two photon 
number states are excited. 

To better probe the nonlinear quantum aspects of this 
dot-cavity coupling regime, in Fig. [3] we increase the 
pump value to r]c = 2.5 g. Here the WEA breaks down 
dramatically, as shown by Fig. [3fa); Figs. [3fc)-(d) fur- 
ther confirm that we are accessing a regime where both 
the WEA and the semiclassical approximations may fail, 
even for a weakly coupled system. We are thus already in 
the anharmonic cavity- QED regime. Figure 3(b) shows 
the corresponding populations. 

Next, we investigate a QD-cavity system in the 
intermediate-to-strong coupling regime, with g = 
60 /ieV ^ tit^ and r]c = 0.25 g. In Fig. HJ^a), we show 
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FIG. 4: (Color online) As in Fig. [21 but for an intermediate- 
to-strongly coupled system, with ^ = 60 /ieV ~ and r]c = 
0.25 g. The other parameters are the same as in Fig. [2l 
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FIG. 5: (Color online) As in Fig. [H but with the larger pump 
field r\c = 1-2 



the transmission with (dashed) and without (solid) the 
WEA; we also show the populations [H^b)], the semiclas- 
sical error and the Fano factor [H^c)]. For this relatively 
weak drive, we again observe noticeable differences in the 
WEA and non-WEA predictions; we also recognize that 
higher-order photon moments are significant already in 
the m = 2 photon state [Fig. 4(b)]. 

With increasing drives, namely when r]c = 1.2 ^, Fig. 5 
confirms that the differences between the WEA and 
multiphoton calculations are even more dramatic (as 
expected), and the on-resonance transmission becomes 
much smaller with the nonlinear drivoi^. Indeed, we rec- 
ognize that multiphoton effects are now fairly profound, 
easily exciting the first 15 photon states. This means 
that we are now accessing the first 30 states of the JC 
ladder, with modest exciton-cavity coupling rates (i.e., 
9 ^ i^t)' 



B. Influence of acoustic phonon coupling 

For the phonon bath calculations, we use material 
parameters for InAs QD^*^, with cj^ = 1 meV and 
ap/(27r)^ = 0.06 ps^. The phonon ME model consid- 
ers a bath at a temperature of T = 20 K, result- 
ing in (5)(20K) = 0.73. We also now consider a 
QD-cavity system in the strong coupling regime, with 
^ = 120 fieV ^ 2 tZf^ using a low cavity pump rate of 
r]c = 0A5g. 

From the theory described in Sec. [Ill electron-phonon 
scattering is seen to manifest in a coherent renormaliza- 
tion in ^ ^ (B) g as well as mediate incoherent scatter- 
ing between the exciton and cavity. To better highlight 
the role of incoherent scattering separately, we first set 
(B) = 1 (i.e., we neglect the coherent phonon effects) for 
the calculations in Fig. [6]^ a) and Fig. [6fb). Figure [6]^ a) 
shows the transmission results with (solid) and without 
(dashed) incoherent phonon scattering. We recognize a 



clear change in the maximum and minimum transmis- 
sion regions and a qualitative reshaping of the spectral 
profile. It is common to try and partially mimic the ef- 
fects of phonon scattering by using an effective 7' to fit 
the data. To highlight the differences with such a sim- 
ple Lorentzian coupling approach, in Fig. [6l^b) we try to 
fit the phonon model by increasing 7' l•67^ which 
naturally broadens the ZPL through an increase of the 
(Lorentzian model) pure dephasing process. However, 
comparing the dashed and solid curves of Fig. [6l^b), we 
see that the features near the center peak and the edges 
of the graph are noticeably different; in particular, a 
broadened ZPL incorrectly widens the Lorentzian tails 
and misses the reduction of the center transmission peak. 
In addition, a change in pure dephasing does not obtain 
the coherent reduction of g {g ^ 0.73g since (B) = 0.73 
at 20 K), which is an important temperature-dependent 
effect that we include in Fig. [6fc); this marked reduc- 
tion in the Rabi splitting will increase with temperature, 
which is important to note especially if trying to fit ex- 
perimental data that is probed via temperature tuning. 
For lower bath temperatures, the on-resonance case can 
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FIG. 6: (Color online) Transmission/reflection characteris- 
tics of a weakly excited strongly-coupled coupled system, 
r]c = 0.15 g, with ^ = 120 /xeV ^ 2 Kf The peak exciton and 
cavity photon populations are both below 0.1 (not shown). 
For these calculations we also consider the effect of acoustic 
phonon scattering at a bath temperature of T = 20 K. (a) 
Transmission with (solid) and without (dashed) incoherent 
photon scattering terms, included through £ph; here we ignore 
coherent renormalization effects, i.e., (B) = 1. (b) Transmis- 
sion as in (a), but now teh non-phonon case has j 1.6 j to 
try and mimic the influence of incoherent phonon scattering, 
(c) Transmission with (solid) and without (dashed) incoher- 
ent scattering and a coherent phonon reduction in ^ ^ (B) g 
with (B) = 0.73. (d) Conditional phase gate, where the larger 
phase changes are (j)^ and the smaller phase change is (j)^ 5 solid 
and dashed lines show results with and without phonon cou- 
pling, respectfully. All other system and material parameters 
are the same as in Fig. [2l and no WEA is made (the semiclas- 
sical approximation breaks down dramatically in this strong 
coupling regime). 
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also be asymmetrioi^. We further note that because we 
are exciting through the cavity mode, the effects of the 
drive on excitation-induced dephasing are significantly 
suppressed, which is in this contrast to driving through 
the QD exciton^i^i^i^. 

Finally, we study a conditional phase gate. We define 
the conditional phase through (p^j^ = (p^j^ — (/>^/^, where 

includes the dot resonance and 0^^^ is the phase 
without the dot. Using a semiconductor micropillar sys- 
tem, conditional phase shifts of around 0.03 rad were 
recently observed^^. Figure [6]^ d) shows the conditional 
phases, with (solid) and without (dashed) phonon cou- 
pling. Near the spectral regions near ±0.1 meV, about 47r 
transmission phase change is possible conditioned upon 
the dot exciton being at the frequency regime — which can 
be tuned for example by applying a field-induced Stark 
shift. In refiection, conditional phase changes of around 
±7r/4 are found to be possible; note that the one photon 
results (i.e., the WEA) tend to overestimate this value. 
Including the coupling to the phonon bath is seen to qual- 
itatively change the phase characteristics, and we stress 
that there is no effective 7' which could mimic the same 
phase trends with phonon bath coupling; the large nega- 
tive phase of transmission actually widens with increas- 
ing 7', instead of narrowing then eventually disappears 
for increasing drives. The phase gate characteristics can 
be optimized further by changing the phonon bath tem- 
perature and by changing the exciton-cavity detuning. 



IV. CONCLUSIONS 

We have presented a semiconductor ME formalism that 
can accurately simulate coherent input /output coupling 
of open-system semiconductor cavity-QED systems such 
as planar photonic crystals and micropillar cavities. We 
investigated the role of quantized multiphoton effects and 
pointed out the possible failure of the WEA, which is 
shown to fail even for low input powers and small mean 
cavity photon numbers (much lower than 0.1). For in- 
creasing field strengths, the possible failure of the semi- 
classical approach is also highlighted. We have further 
shown that coupHng to an acoustic phonon bath causes 
considerable qualitative changes to the light propagation 
characteristics than is modeled by a simple pure dephas- 
ing process. Finally, we used this model to simulate a 
conditional phase gate. 
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